COMBINATORICS OF DISPERSIONLESS INTEGRABLE SYSTEMS AND 
UNIVERSALITY IN RANDOM MATRIX THEORY 

YUJI KODAMA AND VIRGIL U. PIERCE 

Abstract. It is well-known that the partition funetion of the unitary ensembles of random matri- 
ces is given by a r-function of the Toda lattice hierarchy and those of the orthogonal and symplectic 
ensembles are r-functions of the Pfaff lattice hierarchy. In these cases the asymptotic expansions 
of the free energies given by the logarithm of the partition functions lead to the dispersionless (i.e. 
continuous) limits for the Toda and Pfaff lattice hierarchies. There is a universality between all 
^ three ensembles of random matrices, one consequence of which is that the leading orders of the free 

Q energy for large matrices agree. In this paper, this universality, in the case of Gaussian ensembles, 

is explicitly demonstrated by computing the leading orders of the free energies in the expansions. 
We also show that the free energy as the solution of the dispersionless Toda lattice hierarchy gives 
a solution of the dispersionless Pfaff lattice hierarchy, which implies that this universality holds in 
general for the leading orders of the unitary, orthogonal, and symplectic ensembles. 
' We also find an explicit formula for the two point function F„m which represents the number 

of connected ribbon graphs with two vertices of degrees n and m on a sphere. The derivation is 
based on the Faber polynomials defined on the spectral curve of the dispersionless Toda lattice 
^ hierarchy, and ^^-Fnm are the Grunsky coefficients of the Faber polynomials. 
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1. Introduction and Background 



We begin with a brief review of the dispersionless hmits of the Toda and PfafF lattice hiearchies (see 
for examples f21 IMl IS])- Here the Toda hierarchy is the 1-dimensional one with the parameters 
{to, ti,t2, ■ ■ ■}■ These limits are found from the equations for the r-functions of the lattice hierarchies 
by subsituting an asymptotic expansion for the logarithm of the r-functions in the small h limit. We 
also describe the connection of the r-functions of these hierarchies with the partition functions of the 
unitary, orthogonal, and symplectic ensembles of random matrix models with the matrix size N (or 
27V) related to TV = ^ (see for examples [T71[T21[2H])- In particular the asymptotic expansions of the 
free energy, the logarithm of the partition functions, as the size of the matrices approaches infinity 
satisfy the underlying assumption for the asymptotic expansion of the logarithm of the r-functions 
of the Toda and Pfaff lattice hierarchies. Then the leading order terms of the asymptotic expansions 
satisfy the dispersionless Toda and Pfaff lattice hierarchies, which are written in terms of the two- 
point functions denoted by F„„i (see for examples [HlElllH])- There is a rich literature dealing with 
dispersionless integrable systems including the KP and Toda lattice hierarchies, which mainly involve 
studying algebraic and complex analytic aspects of the equations and/or discussing applications to 
2-dimensional topological field theories (see, for examples, [U [211 HSl [Ml [13 ES])- However, 
it seems that there is no paper which connects directly the dispersionless integrable systems with 
combinatorial problems associated with random matrix ensembles. The present paper is to deal with 
this connection for combinatorial problems such as counting ribbon graphs on a compact surface 
of genus zero, and give a universality result for the dispersionless integrable systems which is an 
analogue of that for the random matrix ensembles. 

1.1. The Toda lattice hierarchy. The Toda lattice equation (the first member of the Toda lattice 
hierarchy) is given by 



There exists a. sequence {t^ '. fi ^ 0} of T-functions with tq — 1 which generate the cifi^hfi by the 
formulas 



Here Di is the usual Hirota derivative, i.e. for a variable tk which is a flow-parameter for the k-th 
member of the Toda lattice hierarchy, Dk is defined by 



The Toda lattice is expressed in the Lax representation with a tridiagonal semi-infinite matrix L 




n 



We may then write the Toda lattice equation in the Hirota bilinear form. 
(1-1) ■ = 2r„+ir„_i. 





by 



Bi = [L] 



DISPERSIONLESS INTEGRABLE SYSTEMS AND UNIVERSALITY IN RANDOM MATRIX THEORY 



3 



where [L]_ is the strictly lower triangular part of the matrix L given by 
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The Lax matrix is also written in terms of a shift operator A :— exp(^), that is, with the eigenvector 
(j) = (01, 02, • ■ we have 

{L(j))n = 4>n+l + bn(t>n + An-l^n-l 

= (A + bn + a„_iA"^)0„ = X(j>n. 

The hierarchy of the Toda lattice is defined by 
dL 

— ^[L,Bk], Bk^iL"]^, k= 1,2,3,.... 

Otk 

The T-functions are then the functions of infinite variables, i.e. r(t) with t — (ti,t2, . . .). Then it is 
well-known that each T„-function also satisfies the KP hierarchy [21 16] : 

1 



(1.3) 



hk+i(D)- -DiDk]T„-Tn^O, fc = 3,4,5, 



where D = {Di, \D2, |£^3, • • ■) and /i„(x) with x = D are the elementary symmetric functions of 
X = {xi, X2, . . .) defined by 

(00 \ 00 

E x„^" = E /i„(x)^". 
fc=l / n=0 

In particular, the first equation with fc = 3 gives 

which is the KP equation, that is, the function u logr^^ for each n satisfies 



(1.4) 



dx 



, du d^i 



dt dx'^ ^ dx) ^ dy^ ' 



with X = ti,y = t2 and t = t^. 

It is also well-known that the r„-functions of the Toda lattice hierarchy satisfy the following set 
of equations (see for example [1 ): 



(1.5) 



Dk - /ifc(D) T„+i • T„ = 0, fc = 2, 3, 4, . . . . 



The first equation with fc = 2 gives 
(1.6) 



(Z?2-i^i)T„+i-r„ = 0. 



This equation with (1.1) gives the nonlinear Schrodinger equation (with the change ^2 ^ ^2), i-e. 

with 1^ = ^^-^^ and = . Thus, the nonlinear Schrodinger equation is the second member of 
the Toda lattice hierarchy. 

We now briefly summarize the dispersionless limit of the Toda hierarchy: The key ingredient in 
the dispersionless limit is to introduce a free energy, denoted F, 



(1.7) 



r„(t; h) = exp ( ^F{To, T) + 0{h-^) ] , 
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where /i is a small parameter, and T — (Ti,T2, . . .) represents the slow variables with Tj. = Mj. for 
n > 1. Note in particular that the limit h gives a continuous limit of the lattice structure, that 
is, the lattice spacing has the order 0{h) and the limit introduces the continuous variable Tq, 

hn — > Tq. 

The free energy F{To, T) is calculated from the limit, 

F(To, T) = Imi log [t„ (/j-^T; h)] . 

Then the dispersionless Toda hierarchy can be written in terms of the second derivatives, 

F„-i„ — — , TO, 71 > 0. 

The Fjjin plf^y B.n important role for the dispersionless integrable systems, and they are sometimes 
referred to as the two point Junctions of the corresponding topological field theory (see for examples 

[mniii). 

In the limit ^ ^ 0, the a„, &„ variables in the Toda lattice become 



a„(t) = exp (log T„+i - 2 log t„ + log t„_i) — > e 
d 

^n{t) = 7Tr(l0gTn - logT„_i) > Fqi . 

Utl 



Foo 



Then the Hirota form of the Toda lattice equation (1.1 1 becomes 
(1.8) Fn = e^°° , 



which is called the dispersionless Toda (dToda) equation. Likewise the KP equation (1.4 1 gives 

-4Fi3 + 6F^^ + 3F22 - , 



which is the dispersionless KP (dKP) equation. Also the dispersionless limit of (1.6) is 

Fo2 — 2Fii — Fqi = . 

One should also note that in the dispersioless limit, the spectral problem Lip — Xcj) with the Lax 
operator and the eigenvector (p = {(j)i,(j)2, ■ ■ ■)'^ leads to an algebraic equation (the spectral curve of 
the dToda equation) , 

g^OO 

(1.9) A = p(A)+Foi + — -, 

where p(A) is a quasi-momentum given by limn^Q A(f> / (j) = exp^dS/dTo). This is a spectral curve 
of genus 0, and the dToda hierarchy defines a integrable deformation of the curve (see Lemma 2.1 1. 
Here the eigenvector (p is assumed to be in the WKB form as —> 0, 

0„(t;;i)=expQ5(ro,T) + O(l) 

Although the function S plays a key role for the dispersionless theory, we will not use it in this 
paper. We just mention that the S-function can be expressed in the form (see for example |31jl. 

oo oo 1 

(1.10) S{To, T) = E + To log A - ^^fc(To, T) . 
The dToda hierarchy may be formulated by the form, 

\ -D{\)Fo _ ..p-D{p)Fo 



(1.11) 



^-D{\)D(p)F t g(D(A)+D(^))Fo 



A — /i 

JO 
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Figure 1.1. A (ribbon) graph with one oriented vertex of degree 8 on a sphere 
(left) and a graph with two oriented vertices of degrees 3 and 5 on a sphere (right). 
Examples of complete connected graphs are shown. The edges drawn are to repre- 
sent ribbons laying flat on the sphere. 



where 

oo I g 

This is obtained by the reduction from the 2-dimesional dToda hierarchy in [8j |31] with the con- 
straints tk = ^tk for all k, and the Toda hierarchy with those constraints is referred to as the 
1-dimensional Toda hierarchy. Then one can show that the dToda hierarchy ( |1.11[ ) leads to the 



(1.5) 



curve (1.9) (see Lemma [2. 1 1 below) . We will further discuss the detailed structure of the hierarchy 
(1.11 1 in Section [2] where we show that ( |1.11| indeed gives the dispersionless limits of (1.3) and 



1.2. Unitary ensembles of random matrices. A special class of solutions to the Toda lattice 
equations defines the partition functions of the unitary ensembles of random matrices. Explicitly 
the r-functions are taken to be 

T„(t) = Z(2)(Vb(A); t) = / dAi • • • / dA„ Y[ |A, - A,f exp (- £ Vt{Xk)] , 

R R \ k=l J 

00 

where Vt(Afc) = Vii{\k) — ^^^^ ^(•^) is a polynomial of even degree (e.g. Vo(A) = for 

the Gaussian ensemble). This function represents integration over the eigenvalues of the random 
matrices of the unitary ensemble, and its Taylor expansion in t contains the moments of traces of 
powers of the matrix (see for example a common reference |27^ ) . 

A ribbon graph is a one complex embedded on an oriented genus g surface, such that the com- 
plement of the complex is a disjoint collection of sets homeomorphic to discs. We may think of a 
ribbon graph as a collection of discs (vertices) and ribbons (edges) such that the ribbons are glued 
together preserving the local orientations of the discs in other words the ribbons lay flat on the 
oriented surface. Our ribbon graphs will be labeled in the following way: The vertices are labeled 
so that they are distinct, and the edges emnating from each vertex are labeled as distinct as well. 
We also define the degree of the vertex as the umber of edges attached to the vertex. We represent 
ribbon graphs for this paper as graphs drawn on an oriented surface (we only consider a compact 
surface of genus zero (i.e. a sphere) in this paper, and refer a ribbon graph on the sphere as a ribbon 
graph of genus zero) . In Figure |1.1[ we illustrate one vertex of degree 8 and two vertices of degrees 
of 3 and 5 on a sphere. 

A fundamental result of random matrix theory is that for a quadratic potential Vo(A), log(Z]Y ) 
posseses an asymptotic expansion in even powers of A^, whose terms give generating functions par- 
titioning ribbon graphs by genus: 
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Theorem 1.1 (Bessis-Itzykson-Zubcr [5,, Ercolani- McLaughlin [H]). For the Gaussian unitary 
ensemble (GUE), we have 



log 



^(2) 



N 



A^:iVT 



7(2) 



N 



= E eg(T)iV2-2<, 



where 



es(T) 



0<iij2,... Jl!j2!--- j J! 



T/ie coefficient Kg(j) i^wes i/ie number of the connected ribbon graphs with jk labeled vertices of degree 
k for k = 1,2, . . . on a compact surface of genus g. 



Note that Theorem |1.1| does not contain our Tq = ^ = nh variable in large N limit. The 



Tq- variable can be inserted naturally by taking the r-function in the form with t = A^T (or T = ht), 



r„(t;;.) = Z(2)(^A2;iVT] 



Then the free energy is given by 



log [Tn{h-'T;h)] =log 



We note here that the scaling A VTq^ leads to 



N ^ 



log 



^(2) 



N 



N 



A^ A^T 1 = Zp) 



2To To 



A^-^T] =r,^Zp)r?A2;nt 



where Tj = T^/'^ ^Tj (which is the so called Penner scaling). This expression then has the asymptotic 
expansion including Tq, 

log [r„(;i-iT; h)] = Y: e<,(t)n2-29 + C(To; N) 
g>o 

= E T^-^'cgiT) + dTo; N) , 

g>0 

where C(To;iV) = log[zi^^(f A^; 0)] which can be computed directly (see Appendix A). Note that 
this asymptotic expansion agrees (at leading order) with the assumption (1.7 1; and thus gives a class 
of solutions of the Toda lattice satisfying this assumption. Then the leading order of this asymptotic 
expansion gives a solution of the dispersionless Toda equation, 



F(To,T)= lim -^log 



^(2) 



N 



A^iVT 



T2eo(T) + Co(To), 



where Co(To) is the leading order of C{To; N) for large limit of N, and as shown in Appendix B, it 
is given by 

(1-12) Coin) ^ ^T^ (logTo - 

This formula has also been found directly from the dToda equations in a connection to the topological 
field theory (see (8.6) of p.209 in 4 ). 

In particular, Theorem 1 1 . 1 1 implies that the two point function i^„m(l; 0) for mn ^ represents 
the number of connected ribbon graphs with two vertices of degrees n and m on a sphere (g — 0), 
i.e. 



F„™(1;0) 



-(i;o) 



dT„dT^ 



-(O) = Ko(0,...,l,0,...,l,0,...), 



n,m> 1 . 
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where I's in Ko(j) arc at the n- and m-th places. We will give an explicit formula for i^„„(l;0) in 
Section [3l 

We also attach an enumcrative meaning to derivatives with respect to Tq: 

Corollary 1.1. Each derivative i^o,2fe for k > 1, corresponds to counting the number of connected 
ribbon graphs with a vertex of degree 2k and a marked face on a sphere. 

Proof. We start with our relation 

(1.13) FiTo; T) = T^eo{T) + Co(To) 

j J- 

j J- 

One then notes that the number of vertices of a ribbon graph with vertices of type j — (ji, j2, ■ ■ • ) 
is V — E ji and the number of edges is e — ^ E^ji- If the ribbon graph is of genus 0, then by the 

i i 

Euler characteristic formula, the number of faces is 

i ^ 



SO that by inserting this equation into (1.13), we find 

j j 



(1.14) F(T^- T) = E '«o(j)r(/'4r + ^o(^o) 



from which one sees that differentiating with respect to Tq is related to counting ribbon graphs 
with marked faces in the enumerative part (i.e. except for the contributions coming from Cq). For 
example: -Fo,2fc(l;0) with fc 7^ is the number of ribbon graphs with a single vertex of degree Ik 



and one of the / = 1 + faces marked. These numbers will also be found in Theorem 3.1 
consequence of the dToda hierarchy. □ 

1.3. The Pfaff lattice hierarchy. The Pfaff lattice hierarchy is normally referred to as the DKP 
hierarchy (The "D" stands for the D-type Lie group which is the symmetry group for the hierarchy 
[5n]). The first member of the DKP hierarchy, called the DKP equation, is given by the set of 
equations 



(1.15) 



The left hand side of the first equation is just the KP equation, and the right hand side gives a 
coupling term with the field v^. Because of this, the DKP equation is sometimes called a coupled 
KP equation [18 . 

In terms of the T-functions, u and are expressed as 



(1.16) " = log"^"; V 



^1 ± '^n±l 

dtj 



The DKP equation is then given by 

\ {2D3 + D'fT 3DMrn±i • T„ - . 



(1.17) 
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It is known that the r- functions are given by Pfaffians [TH1[31|52]. Since the equation with the set of 
T-functions has a lattice structure (i.e. t„ is determined by the previous r-functions with tq — I), 
we call it the Pfaff lattice in this paper. This choice of name is also convenient for this paper as 
we wish to compare the Pfaff lattice, as given in [3j |22], and the Toda lattice. The Pfaff lattice 
hierarchy is then given by [21 [31] 



[hk+iCD) - iDiDk+s) T„ • T„ = /ifc(D)T„+i • T„_i, 
^ ^ fc = 0,1,2,.... 

hk+siT'D) Tn±i • r„ = 0, 

In a similar manner as in the Toda lattice case, one can also consider the dispersionless limit of 
the Pfaff lattice equation. The dispersionless Pfaff (dPfaff ) lattice equation is then given by 

(1-18) <^ 

{ 2_Fo3 + Fqi + QFqiFii ~ 3-Foi^02 — 6F12 = . 

Here we have used the limits, = exp(log r„±i— log r„) exp(±/i~^_Fo) E^nd v^v~ = exp(log t„_|_i — 
21ogT„ + logT„_i) exp(Foo) as /i ^ 0. 

In [31] , Takasaki formulates the dispersionless DKP (Pfaff lattice) hierarchy by taking the disper- 
sionless limit in the differential Fay identities of the Pfaff lattice hierarchy: 



(1.19) 



\2 -D(A)Fo _ ,,2 -D(a<)-Fo 
^DWD(,W ^ ^ _ p(^) ^ ^(^))^^ „ ^ 

A ~ /i 



This equation is somewhat similar to the dToda hierarchy ( |1.11[ ), and in fact, this similarity is one 
of our motivations to consider the universality between dPfaff and dToda hierarchies. 

1.4. Orthogonal and symplectic ensembles of random matrices. Two special classes of so- 
lutions of the Pfaff lattice hierarchy are given by the partition functions of the orthogonal and 
symplectic ensembles of random matrices (see for examples [211 [21 [22] ) . We denote those partition 
functions by Z^^'' with P — 1 for the orthogonal ensembles and /3 = 4 for the symplectic ensembles. 
With this notation, the partition function (or r-function) for the Toda hierarchy is denoted by Z^^)^ 
i.e. P = 2. 

1.4.1. Orthogonal ensemble. In this case, the partition function is taken to be 

zW(yo(A);t) = /dAi---/dA2„ ni^'-^^l e^pf- £ ' 

B R V fc=l / 

00 

where Vt{Xk) ~ Vb(Afc) — J2 ^k^j Vb(A) is a polynomial of even degree. This integral represents 

integration over the eigenvalues of the random matrices in the orthogonal ensemble, and as in the 
case of unitary ensemble, its Taylor expansion in t gives the moments of traces of powers of the 
matrix. 

This partition function is also related to the combinatorial problem of counting ribbon graphs 
on a compact surface. In this case we will consider Mobius graphs, which are defined to be one- 
complexes embedded in an unoriented surface of Euler characteristic x such that the complement 
of the complex is a collection of sets homeomorphic to discs. We may think of a Mobius graph as 
a ribbon graph where the ribbons are allowed to be glued together with a twist, i.e. reversing the 
local orientation of the vertices. Our Mobius graphs are labeled in the same manner as the ribbon 
graphs. We introduce Mobius graphs here for completeness sake, however in this paper we consider 
only those Mobius graphs of Euler characteristic 2, such graphs are equivalent, up to a choice of 



local orientations at the vertices, to ribbon graphs embedded in the sphere (see Lemma 5.1), so 
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we will not encumber our discussion here with examples of Mobius graphs. One may see |28j for a 
discussion of Mobius graphs with Euler characteristics less than 2. 

For a quadratic potential Vb(A), we have an asymptotic expansion for log{Z2l^) of the same basic 
character as that of the unitary ensembles: 

Theorem 1.2 (Goulden-Harer-Jackson [ITJ. For the Gaussian orthogonal ensemble (GOE), we 
have 



log 



E 4i)(T)(2iV)> 



X<2 



whe 



TJ 



and /Cx(j) number of connected Mobius graphs with jk labeled vertices of degree k for k = 

1,2,..., on a compact surface of Euler characteristic x- 

We then take our r-function in the form with Tq = — nh, 

r^{h''T;h) = A^iVx) = (2To)"^+tzW (|A^;2nt) , 

where T is given by the Penner scaling, i.e. T ~ (Ti,T2, . . .) with Tj = {2T{))i^^Tj. This expression 
then has the asymptotic expansion 

\og[Tr,{h-'T-h)\ = E 4'nT)(2n)>^ + c«(ro;iv) 

X<2 

= E(2ro)^4^HT)iVX + C(i)(To;iV), 

X<2 

where C^^\Tq\ N) = log[Z2^''(^A^; 0)]. Note that the form of this asymptotic expansion agrees (at 
leading order) with the assumption (1.7). For future reference, for this initial condition, the leading 
order of the expansion gives a solution of the dPfafF hierarchy. 



f(^)(To,T).^lim^^i,log 



j(i) 

-'2n 



N . 



-A^TVT 



= 4T2i?W(t) + cW(To), 

where c''q\Tq) is the leading order of c''q\Tq\ N). In Appendix B, we explicitly calculate Cq^''(To), 
and it gives 



C'o 'in) - To^log(2ro) - -T^ = -Cr(2ro) , 



where Q^(To) = Co(ro) given in ( |l.l2| . The relation with C(2)(To) of the dToda hierarchy is a 
consequence of the universality, which will be discussed in Section [5] 

With Tj = (2To)i~^Tj, we have the expansion formula for the solution of the dPfafF hierarchy 
corresponding to the Gaussian orthogonal ensemble, 

(1.20) F(i)(To,T) = i2TofEi'\T) + C^„'\To) 



TJ 



= (2To)^E^2(j)(2ro)^-''^ + c(^HTo) 



E^2(j)(2To)/^ + C^^)(ro) 



.TJ 
J! 

This formula then gives the combinatorial meaning of the solution of the dPfaff hierarchy (also notice 
the similarity with the formula (1.141 for the dToda hierarchy. 
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1.4.2. Symplectic ensemble. In this case, the partition function is given by 

zi^\Vo{Xy,t) = JdAi • • • /dA„ n - ^jl* exp (- E 214(Afe) ) , 



k=l 



where Vt(Afc) = Vb(Afe) — E ^i^j ^^^d Vo(A) is a polynomial of even degree. This function represents 

integration over the eigenvalues of the random matrices in the symplectic ensemble, and the moments 
of traces of powers of the matrix are obtained in its Taylor expansion in t. 

Then the corresponding combinatorial result is give by the following Theorem: 

Theorem 1.3 (Mulase-Waldron |28|). For the Gaussian symplectic ensemble (GSE), we have 



log 



N 



N 



-A^O 



E 4^Ht)(-27V)> 



X<2 



where e'^\'T) is defined as in Theorem 



1.2 



This Theorem can also be found from Theorem 4.1 of [7J, in particular the 1/(4A^)" appearing 
in that theorem leads to the 2A^T in the above formula of Theorem 11.31 



Remark 1.4. There is a duality in the Gaussian orthogonal and symplectic ensembles. Theorems |1.2| 
and |1.3| give the N ^ —N duality between the partition functions of the Gaussian orthogonal and 
symplectic ensembles of random matrices. In general the even terms of the asymptotic expansions 
agree, while the odd terms have opposite signs. See the discussions in p8j Ff]. 



The r-function in this case is defined by 

'N . 
1 



where T is the Penner scaling defined in the same form as in the case of orthogonal ensemble, i.e. 
t = (fi,f2, . . .) with fj = {2TQ)i-^Tj. Then the free energy is given by 



log [T^{h~^T;h)] =log 
= log 



N 



N . 



-A^2nT 



K4) 



-A^2nT 



+ CW(To;iV) 



= E(-2ro)'^4^nT)iVX + CW(To;iV), 

X<2 



where C^^^Tq; N) — log[Zn\^X'^', 0)]. In the limit N —^ oo, the leading order of the free energy is 
given by 

(1.21) 



i^W(To,T)= hm -^log 

N^oo iV 

^2B^(l) 



(2To)^i?^^^(T) + Cr(To) 



(4)/ 



where Cq (Tq) is the leading ord er of C '^^^Tq; N), and in Appendix B it is found to be the same as 
Cq^\Tq). This implies that with (1.201 we have the duality of the free energies, i.e. 

F«(ro,T)=FW(To,T). 

In Section [5] we further discuss the universality for those free energies corresponding to the GUE, 
GOE and GSE random matrix theories. 
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1.5. Outline of the paper. The paper is organized as follows: In Section[2] wc first show that the 
dToda hierarchy is determined by variables Fqo and Fqi (Proposition |2.1[ ) . Then we discuss some 
properties of F„i„ (Proposition 2.4 1. In particular, we show that Fi„ are related to the Catalan 
numbers, which gives a combinatorial meaning of the dToda hierarchy (Proposition 2.2). In Section 
|3] we first derive the explicit formulae for the two point functions Fnm for the dToda hierarchy under 
the conditions of Fqi = and Fn = 1 (Theorem 3.1 1, and then discuss the combinatorial description 
of Fmn as the enumeration of ribbon graphs on a sphere (Proposition |3.1[ ). We also show that the 
combinatorial meaning of Fmn can be directly obtained from the dToda hierarchy (Theorem 3.2 1. 
In Section [4] we describe the dPfaff hierarchy. In particular, we show that under the assumption 
of the special dependency of the free energy F in the variables (Tq.T), the dPfafF hierarchy can 
be reduced to the dToda hierarchy (Proposition 4.1). This corresponds to the dispersionless limit 
of the PfafF lattice hierarchy restricted to symplectic matrices [23] . In Section [s] we first show the 
universality among the free energies for the GUE, GOE and GSE in the leading order of large N 
expansions (Proposition 5.1). Then we show in general that the solution of the dToda hierarchy 
also satisfies the dPfaff hierarchy under the rescaling of the two point functions Fmn (Theorem 5.2 1. 
This means that the dPfafF hierarchy contains the dToda hierarchy, and shows the universality of 
the dToda hierarchy at the leading orders among the unitary, orthogonal and symplectic ensembles 
in random matrix theory. In Appendix A, we give a brief introduction of the dKP theory, and derive 
the dKP hierarchy (2.3 1 which plays an important role in this paper. We also give in Appendix B 
explicit computations of the terms clf \To) for /3 = 1, 2 and 4. 



2. The dispersionless Toda hierarchy 



In this section, we first show that the dToda hierarchy ( 1.11 1 defines a deformation of the algebraic 
curve given by ( 1.9 ), and as a result we show that the two point functions Fnm are determined by Fqi 
and -Foo- We also discuss a connection of the dToda hierarchy with the dKP hierarchy, and define 
the Faber polynomials as the generators of the higher members of the dKP hierarchy. These results 
will then be used in Section [3] to derive explicit formulae of Fmn solving the two- vertex problem. 

2.1. The dToda curve and its deformation with the dToda hierarchy. Let us first define 
the function p(A), 

p{X) -.^ Xexpi~DiX)FQ). 



This corresponds to the equation p(X) = J^- with the 5-function given by (1.10). Then the dToda 
hierarchy is defined on the dToda curve 

Lemma 2.1. The dToda hierarchy ( |1.11[ ) gives a deformation of the algebraic curve of genus given 
by 

X^P + Foi + . 

P 

Proof. Eliminating e^'^^^^'^^^^ from two equations in ( |1.11[ ), we obtain 

V A-/i )\ p{X)p{t,))- 

We then separate the variables, 

A -_p(A) T-- = p-pifj.) 



pW PifJ-)' 
With the asymptotic condition, A — p(A) Fqi + 0{X~^), we obtain the curve. □ 
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With this Lemma 2.1 we define in this paper the dispersionless Toda (dToda) hierarchy in the 
form, 



(2.1) 



Foo 



A = p(A)+Foi + — with p(A) = Ae^^(^)-^° . 

pW 



In particular one should note through the paper that the dToda curve (the second equation) on 
(A,p)-plane plays a central role for the combinatorial description of the dToda hierarchy. 
We now note that the first equation in (2.1 ) implies: 

Lemma 2.2. Each Fmn with m,n > I can be determined by the set {i^ofe : < A: < m + n}. 

Proof. The first equation with p{X) = A exp(— £'(A)^o) gives 

The right hand side depends only on Fq^ for A; > 0. □ 

Also, setting A = /i and expanding for large A, the first equation gives 

pFoQ pFoQ oo 1 

where Fq = (-Fbi, ^02, • • •) with Fok ~ ^Fok, and Z = (0, Z2, Z3, . . .) with Z„ for n>2 defined by 

Zn= E ^ (A:,/>1). 
k+l=n '>^' 

Comparing the coefficients of the power A~", we obtain 

(2.2) /i„+2(-Z) + e^»/i„(2Fo) = n > 0. 

This set of equations is called the dipersionless Hirota equation for the dToda hierarchy, which gives 
a half set of the dispersionless Hirota equations (see for example |8] , and the other half is given by 
(2.4) below). The first equation with n = then gives the dispersionless Toda equation. 

For n = 1 , we obtain 

i^i2 = 2e^™Foi=2FoiFii. 

Thus Fii and F12 are determined by only Fqq and Fqi. In fact, we have from the second equation 
(dToda curve) that 

Proposition 2.1. The two point functions Fnm of the dToda hierarchy is determined by Fqq and 
Foi only. 



Proof. The curve gives 
Expanding this with large A, we have 



1 _ „-D(\)Fo , , ^pDiX)Fo 

^ A ^ A2 ^ 



/i„(-Fo) + e^''«/i„_2(Fo) =0 n>2. 
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Then one can see that F^n is determined by the previous i^ofc for < fc < n. For examples, we have 
for n = 2 and 3, 

Fo2 - i^oi + 2e^™ = + 2Fn , 

Fo3 = ^FoiFn2 - ^F^i + SFoie^"" = F^, + GFoiFn . 



Then from Lemma 2.2 we conclude that all F„in for m, n > 1 are determined by Fqq and Fqi (or 
Foi and Fn = expFoo). □ 

From the dToda curve, p{X) can be explicitly calculated as follows: 

Proposition 2.2. The inverse of the dToda curve for the case p{X) — > A as A ^ cxd is give by 
where 

Note here that Fn = e^"" and Ck is the k-th Catalan number with Cq — 1. 
Proof. From the curve, first we have 

p2-(A-Foi)p + Fii-0. 
Solving this for p with the asymptotic condition p{\) — > A as A — > oo, we have 

p=\{\-Foi + v/(A - Foi)2 - 4Fn 



2 

= A - i^oi - ^ (a - ^^01 - v/(A-i^oi)2-4Fu) 
Noting the last term can be expressed as 

Fii 1 - yr^4^ Fn ^ ^ k ^^11 



J2 CkX with X = 



A-i^oi 2a: A - Fqi (A-^^oi)^ 

Thus we have 



Then the coefficients P„ can be found by 



Fi,„+i _ - ^ ^,+1 r dX (A + foi)" 



Using the expansion (A + Fqi)" — ^ (^)X^ F^^ , we obtain the formula. □ 

The appearance of the Catalan numbers indicates some connections of combinatorial problems to 
the dToda hierarchy. This is one of the main motivations of the present study. 

2.2. The dToda hierarchy implies the dKP hierarchy. Here we consider the second equation 
of that is, in terms of p{\), 

(2.3) pW-pM ^^d(x)dmf_ 

A /i 

This equation is also known as the dKP hierarchy, and p{X) in the dKP hierarchy is given by 
p(A) = A — D{X)Fi (see Appendix A, and also for examples [S] |S]). For the dToda hierarchy, 
p(A) = Aexp(— D(A)i^o) can also be expressed in a similar form. This is obtained from an interesting 
connection with the Faber polynomials: 
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Definition 2.1. Let A be a Laurent series in p given by 

X^p + ui + }^ — ^. 

fc=i P 

Then the Faber polynomials ^n{p) of degree n in p are defined by 

^n{p) ^ [X{pr]>o n = 0,l,2,..., 

where [A(p)"]>o represents the polynomial part of A(p)" in p. The polynomial <i>„(p(A)) has the 
Laurent series in A, 

oo 1 

<i>„(p(A)) = A" - E —Qnm. 

m=l 

where Qnm are the functions of u (ui, U2, • • ■)■ Those Qnm are related to the Grunsky coefficients 
Cnm = -^Qnm- Then (2.3 1 gives the generating function of the Faber polynomials with Cnm = j^Fnm 
(see Remark 2.2 below for a further discussion). In the dKP theory, the Faber polynomials give the 
generators of the flows in the hierarchy (see Appendix A) . 



Then with (2.3 1, we have 
Proposition 2.3. The Faber polynomials can be expressed as 

<i>„(p(A)) = A"-i?(A)F„ n>l, 

and ^o{p) = 1- 

Proof. We first note that 

( pw^pii^y 



log [^^Yzf^ ) = Di\)D{^)F. 



Taking the derivative with respect to A, we have 



p{\)-p{^i) d\ X-^i ,^iA»+ 
Then for n > 1, we have 

/ dX A" dp(X) „ , , 
/p=oo 27ri p - p{^i) 

This completes a proof. □ 

We now derive the other set of the dispersionless Hirota equations: Taking the limit X ^ p oi 
the equation in (2.3 1, we have 



dp{x) 



dX 

With p{X) = Aexp(-i:)(A)Fo), we have 

-X^D{X)Fo = e^(^)'^+^(^)^o _ 1 



dX 



D 1 \ OO 1 

\n=l J n=l ^ 



where Z = (Zi, Z2, . . .) with Zk defined by 

Zi-i^oi, Z„-— + 2^ -rr. 

k+l=n 

This then gives the other half set of the dispersionless Hirota equation for the dToda hierarchy (the 
first half set is given by ( |2.2[ )), 

(2.4) Fon^K{^) n>2, 
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which gives the dispersionless limit of the hierarchy (1.51. The first equation of (2.4 1 with n = 2 
gives 

^02 = -Fqi + 2i^ii . 

With Proposition 2.3 with the case n = 1, i.e. ^i{p) — p + -Foi, the function p{X) can also be 
expressed by 

p{X)^X-Foi-DiX)Fi. 
Then note that that the dHirota equation of the dKP hierarchy is derived from the same equation 
dp/dX = e^(^)'-^ with p{X) = A - Fqi - D{X)F, 

Fin = K+i{Z) n > 3, 

which is the dispersionless limit of ( 1.3 1. Here recall that Z = {0, Z2, Z^, . . .) with Z/^ is given by 



Zk = Zk — Fok/k. The first equation with n ~ 3 gives the dKP equation 
(2.5) - 4Fi3 + 6Fi\ + 3F22 = . 

We also mention the positivity of Fmn'- 

Proposition 2.4. The coefficients Fmn satisfy the following properties: 

(a) // Fim > for all m > I, then Fmn > for all m,n > 1. 

(b) If Fi^2k = for all k > I and all others i^i.m > 0, then Fmn = for all m, n with m+n =odd 
Proof With p{X) ^ X- Fqi- D{X)Fi, we have 

^ p{X) ^ pjfi) ^ ^ {D{X) - D{^i))F 
A — /i A — /i 

1 00 r / 1 1 



1 00 P / m 1 



Xfi m=i m Vfe^i A™ fc^*^ \ 
Then taking the logarithm and the differentiation with respect to A of this equation, we have 

Fmn 1 ^ -D(\)Dit.)F ^ Fim f ^ m - k + 1 



m n 



This shows the property (a), i.e. VFim > implies ^Fmn > 0. 

The second property (b) can be obtained from dp{X)/dX ~ e^^"^^ ^ . This equation gives 



k 



=1 Afe^ 



which shows that if all Fi_2k = 0, then 



n+m— odd 



Then from (a), we conclude Fmn = if m + n =odd. □ 

Remark 2.2. The Faber polynomials <&„(p) appear in a classical complex function theory: The 
Bieberbach conjecture sates that if f{z) is a univalent function on \z\ < 1 with the expansion 
f(z) = z + a2Z^ + a^z^ + • • • , then the coefficients a„ satisfy |a„| < n. Let g{z) be the function 
defined by 
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which is univalent for \z\ > 1. Then the Faber polynomials <&„(p) appear in a expansion, 

\ z J nz" 
With the expansion of <i>„((7(z)) for large z, 



(2.6) ^r^{g{z))^Z^ + 7lY. 



m—1 



z'- 



Grunsky defined the Grunsky coefficients c„„i- Then Grunsky showed that his coefficients should 
satisfy a sequence of inequalities in order that g is univalent on \z\ > 1. In particular, an explicit 
formula for the Gransky coefficient was found in the paper by Schur in terms of the coefficients a„ of / 



[29] . In the section 3.1 we find an explicit formula of Fnm for the case where the cefficients are given 
by the Catalan numbers and discuss their combinatorial meaning (our derivation is independent 
from the calculation used in [29j). 



From (2.3) and (2.6), one can see that the Grunsky coefficients Cnm are related to our two point 
functions, that is, Cnm = j^Fnm- The dToda hierarchy then defines an integrable deformation 
with an infinite number of parameters T = (Ti,T2, . . .) for the univalent function f{z), so that the 
Grunsky coefficients are given by the second derivatives of a function F, the free energy, i.e. 

1 d^F 

nm OlnOlm 

Thus a most important consequence of the dToda hierarchy is to show the existence of a unique 
function F , the free energy. 



3. COMBINTORIAL RESULTS FOR THE DTODA HIERARCHY 



The appearance of the Catalan numbers in the dToda hierarchy (see Proposition 2.2) directly 



demonstrates a combinatorial meaning of the functions Fnm- In particular, if Fqi{1\0) — and 



Fii(l; 0) — 1, then we have, from Proposition 2.2 

P(A) 

which gives 



(3.1) 



Fi,2fc+i = (2fc+l)Cfc 



In general, the two point functions Fnm have the following combinatorial interpretation: One sees 
from the connection to random matrices discussed in section |1.2| that for the particular choice of 
F(To, T) = ro2ep(t) + Cn(ro) the two point function with n, m > 1, 



^^nm(To = l;T = 0) = 



d'^F 



dTndT„ 



92 



= 1:T=0 



dTndTr. 



-7]?eo(T) 



= 1:T=0 



is the number of genus connected ribbon graphs with a vertex of degree n and a vertex of degree m 
(see Theorem 1 1.1| . Here the degree of the vertex represents the number of edges (ribbons) attached 
to the vertex. For example, the i^i,2fc+i gives the number of connected graphs with a vertex with a 
single edge and a vertex with 2fc+ 1 edges on a sphere (recall that the number of ways to connect 2k 
edges of a single vertex is given by the Catalan number Ck) ■ In this section we compute a remarkable 
closed formula for the i^„,m(l; 0) in the case with ^oi(l; 0) = and -^11(1; 0) = 1. Also recall that 
differentiation with respect to Tq has the meaning of counting ribbon graphs with a marked face. 
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3.1. Explicit formulae for i^„m(l;0). Let us first recall Proposition 2.4 (b): if Fi 2k — 0,Vfc > 
and i^i.m > OjVto > 0, we have 



F — n 

71,171. ^ 



if n + m = odd . 



Also from Proosition 2.2 if Fio = and Fn = 1, we have Fi^2k = and Fi_2fc+i = (2fc + l)Ck- So 
we calculate only F„m with m + n =even and positive. Then we state the following Theorem which 
gives an explicit formula of Fnm for those cases of m and n. The method to find the formula is 
based on the dispersionless Toda hierarchy. 



Theorem 3.1. The two point functions F^m — qt qt ^ f'^^ dToda hierarchy (2.1) with J^oi = 
and Fii = 1 (i.e. Foo = are given by 

Fo,2k = (fc+l)Cfe, fc=l,2,..., 



F 



(2j + l)(2fc + l)^-^^^>^C,Ck, J, fc = 0, 1, 2, . . . , 



2j + l,2k+l 



Fo 



.(J+l)(fc+l) 



j,fc = l,2,. 



Fnm — 0, otherwise , 
where Cj is the j-th Catalan number. 

Proof. We first derive the formula Fo,2fc for k ^ 0: We start with the definition of p, 

p(A) = Ae-^(^)^° . 

which implies 



Therefore we may compute by the contour integral 



1 

2k 



Foak — 



A=oo 27ri 



A^^-^log 



P(A) 



Using the dToda curve, 



\{p)^P+ - =p\ 1 + 



P 



we write the integral with respect to p, and then integrate by parts, 



1 

2k 



Foak — 



dp 1 
p=oo 27rj 2k \dp 

p=oo ^TTJ k 



\{pY^ ) log ( 1 + 



1 



2fe-l 



2k-l 



P" 

I (2k -I 
k\k-l 



In the last step, we have used the expansion, (1 + ^ = ■ Thus we obtain 

J^0,2fc = (fc+l)Cfc. 



i.e 



We now derive the formula Fnm for nm ^ 0: We use the Faber polynomials in Proposition 2.3 

a>„(p) = [A(p)"]>o = A"-i^(A)F„. 

Fmn f dX 



This implies 



m 



27ri 



A'»-^<i>„(p(A)) . 
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Now changing variables from A to p, we have 



dp 
27ri \^ dp 



dp 
2Tii 



^ Pj dp 



P+ - 
P 



Here we have noticed that ^n{p) can be replaced by A(p)" with (i([A(p)™]<o]) for d(A(p)™). The 
notation [/(p)]<o implies the part of the Laurent series /(p) with the inverse powers of p, i.e. 
[/(p)]<o = f{p) - [f{p)]>o- Thus we have 



(3.2) 



= I] (to - 2fc) 



fc=0 



dp 

2TTi 



J2k-m-l 



It 

Then from the expansion [p + i)" = ^ we have the residue integral 



dp ^n-m+2{k-l)-l 

which implies n — m + 2(A: — Z) = Oin the non-zero case. Therefore, Fmri 7^ only if n — m is even 
(recall Proposition 2.4 1. So we have two cases: (a) to = 2j + 1, = 2fc + 1, and (b) m = 2j, n = 2k. 
With the residue calculation for the case (a), we obtain 



F2j + l:2k+l = X] (2^ + 1) 
i=0 



2j + l\ /2A: + 1 



k-l 



We consider the case (a) (the case (b) can be shown in the similar way): Let us assume k > j 
(because of the symmetry in j, k, this case is enough to show). Also confirm that the case of j = 
gives (3.11. Then from ( |3.2| , we have 

2j + l\f 2k + 1 

k-j + l 



F2j+i,2k+i - E(2j-2Z + 1) 



We note the following identity which is a key equation in the proof, 



(3.3) 

Also we write 



(2j-2/+l) 



2j + 1 



2j 



2k + I 



= (2j + 1) 
k{k~l)---{k-j + l+l) 



2k+l\k-j + lJ {k + 2){k + 3) ■ ■ ■ {k + j - I + 1) 
Then we have the following formula with the common denominator (fc + 2)(fc + 3)---(fc + j + l). 



F 



2j + l,2fc+l 



3 

(2j + l)(2fc+l) "%?o 



2j 
I - 1 

Cfc 



(fc + 2)(fc + 3)---(fc + j + l) 
where G(j, k) is the polynomial of degree j in k 



k{k-l)---{k-j + l+l) 
(fc + 2)(fc + 3)---(fc + j-/+l) 

G(j, k) , 



Ck 



(3.4) 



G{j,k) :- E 

i=0 



2j 
l-l 



a=0 



i-1 



a X 



JJ(fc + j + 1 -a) 



a!=0 



DISPERSIONLESS INTEGRABLE SYSTEMS AND UNIVERSALITY IN RANDOM MATRIX THEORY 19 

-1 



Here the products JJ^ (• • • ) are understood to give 1 when I = and I = j. We then claim that the 

a=0 

polynomial G{j, k) can be expressed by the simple form, 



(3.5) 



Here note that the coefficient ( r*) is just J2 ( /) 



(fc + l)(fc + 2)...(fc + j) 

(j + l)Q(fc + l)(fc + 2)...(fc + j) 

j 



1=0 



This then gives the desired formula, 



that is, this equation for G{j, k) leads to the formula F2j+i^2k+i- 

To prove the claim, we write (3.4 1 in the sum of the terms with the coefficients (^^^ (like a 
telescoping sum): 



/2j^ 



1=0 



(3.6) G(.7,fc) = -2EJj-0(7) n (fc-a)xn(fc+J + l-«)+(^^)(fc+2)(A;+3)...(fc+j + l), 

a=0 Q=0 ^ 



where the last term in (3.6 1 is the remainder from the telescoping. We then observe that this reduces 

Dec 



the degree of the polynomials of both sides of (3.5 1, that is (3.5) becomes the equation 



1=0 



(3.7) 2E(j-0( ;) n (fc-«)x n(fc + J + l-")=J )(fc + 2)(fc + 3)---(fc + j), 

Q=0 a=o ^ 



where the right hand side of this equation is the difference between the remainder in (3.6 1 and the 
right hand side of (3.5 1. We then note that the coefficients in the sum can be written in the form 
similar to (3.3 1, 



= J 



2j - 1\ f2j - 1 



/ - 1 



Here we have used the relations (the second one is the Pascal rule). 



2j 



= 2j 



2j-l 
/ - 1 



and 



2A f2j - 1 



2j-l 
/ - 1 



Then we can write (3.7) in a similar form to (3.5), 

]-l-2 



i-1 



E 



2.7 - 1\ /2j - 1 



I 



l-l 



i-i 



2j-l 
J-1 



(fc + 2)(A: + 3)---(/c + j) 



II {k-a) X l[{k + j + l-a) 

a=0 a=0 

Again we write the sum in the terms with the coefficients (^"'^^) and move the remainder to the 
right hand side to arrive at the equation 

j-l-3 l-l _ 

' \k + 3){k + A)---{k + j). 



'j:{2j -21-1) 

1=0 



2j-l 



J] {k-a) X l[{k + j + l-a) = 

a=0 a=0 

Both sides are now polynomials of degree j — 2 in k. In the similar way with the identity, 



2j-l 
J-1 



(2j -21-1) 



2j-l 
I 



(2j-l) 



2j - 2\ f2j - 2 
I - 1 



the above equation takes the similar form as before, i.e. 



i-2 

E 

1=0 



2j - 2\ f2j - 2 



I 



l-l 



3-1-3 



l-l 



Yl (k-a) X l[{k + j + l-a) = 



2j - 2 
J-2 



(fc + 3)(fc + 4)---(fc + j). 

\ / — X, / 

a=a a=0 

Then we can further reduce the degree of the polynomials in the same way. Continuing this process, 
we complete the proof. □ 
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2s 



2s+i ; 





2j-2s-l 

Figure 3.1. Two- vertex problem with degrees 2j + 1 and 2k + \. The gluing of 
the edge 1 of 14 to the edge 2s + 2 of the same vertex divides into the two sets of 
edges; the set of 2s edges marked by {2, 3, . . . , 2s — 1} is of the one- vertex problem, 
since those edges can not cross the glued ribbon between the edges 1 and 2 s + 2 on 
a sphere. The gluing of the 2j — 2s — 1 edges of the vertex Va to 2fc + 1 edges of Vf, 
gives the two- vertex problem. 

3.2. Graph enumeration on the sphere. We recall that Fnm has a combinatorial meaning of 
enumeration of connected (ribbon) graphs on a sphere. The presentation of such combinatoric 
problems in terms of the function F{Tq\ T) arising from the Gaussian unitary ensemble example is 
not a new idea. However viewing the genus zero solutions as arising purely from the dToda hierarchy 
is a new approach to solving the problem, although a similar idea was used in the case of a single 
even time for the computations of |16j . Surprisingly we have found a closed form for the numbers 
Fnm with the conditions _Foi = and Fn — 1 for general n and m. 

The solution of the two- vertex problem given by Theorem |3.1| is new, previously the problem has 
been solved only in the so called "dipole" case when the vertices are of the same degree (i.e. F„„) 
[19l|25]. These authors consider the general genus problem of two vertices of the same degree on a 
genus g surface. Similar formulae were also found for other low genus problems in [16|,l5|. but again 
the vertices must all have the same degree. 

that is, 



Let us start by giving a combinatoric argument for the F^ ^k formula of Theorem 3.1 



-Fo,2fc gives the number of genus ribbon graphs with a vertex of degree 2k and a marked face: 
As shown above counting graphs with marked faces in the case of fixed genus is equivalent to just 
counting the number of graphs, the number of faces being given by the vertex structure. So the 
problem is then to count the number of genus ribbon graphs with a single vertex of degree 2k. It 
is easy to see that this is equivalent to counting the number of ways to make pair-wise connections 
with 2k points on a circle without crossing and the number of sections divided by the chords giving 
the connections of the vertices. The first number is well-knwon and is given by the fc-th Catalan 
number Cfc. Then the sphere is divided up to fc + 1 sections. So we have (fc -I- l)Ck which is fo,2fe 



(see Theorem 3.11 



Now we consider the two vertices problem, and show that F^m satisfy the following recursion 
relations as a direct consequence of the enumeration of connected ribbon graphs appearing in the 
problem: 

Proposition 3.1. The two point functions F^n with nm ^ for the two-vertex problem satisfy the 
recursion relations, 



F2j+i,2k+i = 2 E Q_/i^2/-i,2fc+i + {2k + l)Cj+k j, k = 0,1,2, 



1=1 



F2j.2k — 2 ^ Cj^lF2l^2,2k + 2kCj + k~l j, k — 1,2, 

1=2 



Proof. We consider the case of i^2j+i,2fe+i (the case of F2j.2k can be shown in the same way). 
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We will use Va and Vi, to refer to the vertex of degree 2j + 1 and 2fc + 1 respectively. We start 
with the graph, in which a specific edge, say 1, of Va is glued to an edge of the vertex Vf,. The 
number of ways to accomplish this gluing is 2fc + 1 (the number of edges at Vf,)- Having glued the 
edge 1 of Va with an edge of Vt, we have the situation of one- vertex problem of degree 2j + 2k. So 
the number of connected graphs for each situation is given by the Catalan number Cj+k- Thus the 
total number of the edge 1 gluing is (2fc + l)Cj+fc, which is the last term in the equation. Now we 
calculate the other cases where the edge 1 is glued to another edge of Va- Let the edge 1 be glued 
with the edge 2s + 2 (s = 0, 1, . . . , j — 1). This divides the edges of Va into two sets one of which is 



cut off the rest of the graph (see Figure 3.11. This group of edges is marked by {2, 3, . . . , 2s + a}, 
and these 2s edges of Va must be glued together within this group. This is the one-vertex problem 
of degree 2s, and the number of gluings is Cg. Now we are left with the two vertices problem with 
the vertex of degree 2j — 2s — 1 for Va and the vertex of degree 2fc -I- 1. The number of connected 
graphs of this problem is, of course, i^2j-2s-i,2fe+i- This situation with one-vertex problem of degree 
2s and two-vertex problem of degrees 2j — 2s — 1 and 2k + \ appears also the case where the edge 
1 is glued with the edge 2] — 2s -I- 1. Here the group of 2s edges in the one-vertex problem is 
{2] — 2s + 2, 2j — 2s + 3, . . . , 2] + 1}. So we have the total number of connected graphs is 

F2j+l,2k+l = 2 X] CsF2j~2s-l.2k + l + (2fc + l)Cj+/j . 

Then changing the variables in the sum with s = j — I, we have the equation. 

We also note that these relations in Proposition |3.1| define a recursion relation of the Faber 
polynomials on the dToda curve: We consider only the first equation (the second equation can be 
shown in the same manner). Multiplying (2fc+i)A^'°+i ^^"^ ^'"^^ equation in Proposition 
summing up for fc = 0, 1, 2 . . ., we obtain, after rearranging the terms, 

(3.8) <i>2,+i - 2 t C,^i^2i-i = A^^*! - t . 

1=1 1=1 

oo oo 

(Recall that <&„ = A" — J2 TriJ^^™^ ^^^^ $1 = p = A — J2 x^^+i Ck-) This can be considered as 

va=l k=0 

a recursion relation of the Faber polynomials defined on the dToda curve A = p -I- ^, that is, with 



3.1 



and 



(3.8), one can define <&2j+i from $1 = p. Conversely, starting from this recursion, one can derive the 



first equation of Proposition 3.1 This implies that we have an alternative proof of this Proposition 
using only the Faber polynomials, that is, the explicit form of two point functions Fmn obtained in 
Theorem |3.1| satisfy the recursion relations of Proposition |3.1[ 



Let us now show that the relation (3.8 1 can be directly derived from the dToda curve. We first 
note that the left hand side is a polynomial in p. It is indeed easy to see that the projection of the 
right hand side on the polynomial part gives the polynomial in the right hand side. This means that 
the non-polynomial part in the right hand side must vanish, when (A,p) is on the dToda curve, i.e. 
X = p + ^. So we want to prove 



(3.9) [m'^p] 



t c.^iHpr-' 
1=1 



where [/(p)]- is the non-polynomial part of the Laurent series f{p). With the dToda curve, this 
leads to 



fc=o \ / ;=o k=Q V 



Then from the coefficients of the power p '^^^^ and using the Pascal relation, (^^-') — {^■'^ ^) + 



^^\^), we have the following Lemma which proves equation (3.9 1 
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Lemma 3.1. For j > k, we have 

k-ij V k-i 

Proof. Let us first note that we may write the second binomial coefficient as 

2(j -I) -I \f2n-V 

where we denoted n = j — I and i ~ j — k. We also note 
(3 '^^"^ 1\ _{n + i) i f 2n 



2i n\n — i 

Now let z{t) be a generating function of the Catalan numbers, 

oo 

i=o 

which satisfies 1 = z — iz^ with z = 1 at f = 0. Then one can show that 

n=i n\n-ij 

This result can be found by computing the contour integral, using t ~ and dt = ^^dz 



dt 1 



^0 2m 



(z-ir 



dz 



^1 2m (z - 
du 



(2-z) 



2n- 1 

n — i 



i ( 2n 



2n~ 1 



n — i — \j n\n — i 



Then from (3.101, we see that 

^2n- 1 



— — \f(z-iy] 
t=o 2m 2idt ^ ^ ^ ^ 

— ^—(z-iy-'iz-l + tz') 
t^o 2m 2i"+i ^ ' ^ ^ > 

1 (z - 1)* 
t^o 27ri 2-z 



This implies that 



2-z 



Also with Co = 1, we have z ~ \ — 'Y^ Cit^, so that we see by setting i — j — k, 

1=1 



f^—iV °° °° / 2?? — 1 



tn+l 



2-Z i=in=j-k \n-{j-k)^ 

2{m - /) - 1 



oc rn—j+k 
m—j — k-\-l l — l 



(m-l) - {j - k) 
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Then the sum in the Lemma is given by the coefhcient oi P , i.e. 

dt 1 , ^ {z- ly-'' 



- --t(z - 1)- 



27ri (z - 
k-1 



2- z 



u=0 



□ 



This Lemma proves the relation (3.9 1 which leads to the recursion relation of the Faber polyno- 



mials (3.8 1 defined on the dToda curve. The point here is that the combnatorial meaning of the 

given in Proposition |3.1| can be directly obtained from the structure of the 



two point functions F„ 

Faber polynomials defined on the dToda curve (i.e. independent from the random matrix description 
of the two point functions). We summarize this as the following Theorem: 



Theorem 3.2. The Grunsky coefficients c„ 



-Fnm of the Faber polynomials defined on the 

that is, they have 



3.1 



dToda curve A = p + - satisfy the recursion relations given in Proposition 
the combinatorial meaning of counting the connected ribbon graphs for the two-vertex problem on a 
sphere. 

4. The dispersionless Pfaff hierarchy 

In this section, we consider some basic properties of the dPfaff hierarchy. In particular, we show 
that the dPfafF hierarchy can be reduced to the dToda hierarchy by choosing a particular section of 
the solution space which depends only on the even times of the hierarchy (i.e. the symplectic Pfaff 
lattice hierarchy studied in [23] ) . We also discuss a rescaling of the variables based on the structure 
of the r- functions of the Toda and Pfaff lattice hierarchies. That structure is that one is given by 
determinants and other by Pfafhans. 

4.1. The dispersionless limits of the symplectic Pfaff lattice hierarchy. Let us first note 



that the dPfaff hierarchy ( 1.19 ) reduces to the dToda hierarchy as a special choice of the variables 

Proposition 4.1. Let a = \^ and P = fi'^. Suppose F2„i+i 
s„ = 2t2„ for n > I and Sq 



aJ^^ — — for all m > 0, and let 

dt2m+l ■' — ' 



tQ. Then the dPfaff hierarchy (1.19) is reduced to the dToda hierarchy 



^D(a)D{p)F _ 



^-D{a)D{p)F 



a/3 



where F{sq, si, • ■ •) = F(io, ii, • • •) and D{a) = ^ 



-(5 

d 



-omFo 



fc=l 



Proof. The proof is based on the following observation: With a — and -F'2m+i — (Vm > 0), 

1 dF 



D{X)F^ E 



1 2mA2™ dt2m 
_ g 1 dP _ 
ma'" dsjn 



b{a)F. 



It is then immediate to see that (1.191 with i^i = gives ( |1.11| ). □ 

This can be considered as the dispersionless limit of the Pfaff lattice restricted to symplectic 
matrices (a continuous version of the SR algorithm [12]). The dPfaff equation is then given by 



F22 = 46^"" which is Fn 
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The second equation in the dPfafF (dDKP) equation (1.18 1 becomes trivial. 

4.2. The dPfaff curve and the rescaled dPfaff hierarchy. We consider the following scalings: 
(4.1) To^2To, T^2T, and F 



1 

2^- 



The scaling in the time variables may be considered as the size of the matrix models, i.e. the 
matrices we consider for the orthogonal ensembles are twice the size of matrices for the unitary 
ensembles. The scaling in the free energy function is due to the structure of the r-functions, i.e. the 
determinant for the Toda lattice and the Pfaffian for the PfafF lattices (a Pfaffian is a square root 
of a determinant). We will show in Section [s] that those scalings appear in the particular choices of 
our T-functions corresponding to the Gaussian ensembles discussed in the Introduction. 

With the scalings (4.11, F„j„ is scaled by 2i^,„„ for all m,n > 0, and we redefine the dPfaff 
hierarchy (1.191 in the foUwoing form similar to the dToda hierarchy, 

„2Foo 



(4.2) 



< 



1 



g(A)2=p(A)2H 



^ ^-2D{X)D{t,)F 



q{\) ~ qifi) 



X 



Go 



with p{X) = Ae-°(-^)^°, g(A) = A - Fqi - 2D{X)Fi, and 

G2 — F02 — Fqi — 4Fii. 

We call the algebraic equation of {p,q) in (4.2 1 the dPfaff curve. The first equation is the scaled 
version of the second one in (1.19), and the first one in (1.19) is derived from this definition (4.2): 



Lemma 4.1. From (4.2 1, we have 

p{Xf - 



X 



= e^^W^(^)^(g(A) + ,(M)), 



which is the first equation in (1.19) with the scaling Fm, 
Proof. From the dPfaff curve, we have 

g(A)2-9(/i)^=(p(A)2-p(^)2) (i- 



IF 



= (p(A)2-p(^)-,e 
which leads to the equation in the Lemma. □ 



p{x)Mf^)^ 

-2D(\)D{t,)F ( - lM 

X — II 



This shows the equivalence between (1.19) and (4.2). 



Let us now discuss the functional relations among F, 
for the dPfaff hierarchy: We first note 

Lemma 4.2. Each Fmn for m,n > I can he determined by the set {-Fofc 

Proof. From the dDKP curve, q^ 
be determined by {Fq^ : < k < m 



i.e. the dispersionless Hirota equations 



P 



2DiX)Diii)F^-log(l 



< fc < m + n}U {Fii}. 

G2 + e^^oo/p^, we can see that each Fi „i for m > 2 can 
1} U {Fii}. Then from the first equation, we have 

^ / iD{X)-Diii))F, 



X — fi 



Pix)M,)\' ' '^^^'^ 

This implies the assertion. □ 

In the similar way as in the dToda hierarchy, the dispersionless Hirota equations can be derived 
as follows: From the dPfaff curve, we have for n> 1 

/i„+2(-2Fo) + e2^°°/jn-2(2Fo) 
(4.3) 4 „ 4 



1 



Fi,i 



+1 



-FoiFin 



.4 y: 

j+k=n i k 



DISPERSIONLESS INTEGRABLE SYSTEMS AND UNIVERSALITY IN RANDOM MATRIX THEORY 25 



where Fq — {Fqi, Fq2, ■ ■ ■) with Fqj = -j^. This gives the first half of the dispersionless Hirota 
equations for the dPfaff hierarchy, which shows that Fi„ are determined by {Fom : < m < n+1}. 
The first nontrivial equation with n 1 is given by 

1 23 



This is the second equation of the dPfaff lattice equation (1.18 1 with the scaling F„,„ — > 2F„„j. 

The second half of the dispersionless Hirota equations is given by taking A = /i in the first 
equation, i.e. 

_ e^^»» ^ -2D(XfF<k 

Expanding for large A, this gives 

p2i^oo 00 1 

OQ I - \ / 00 

E )^/^«(-2Z)) 1 + 2 E ^^^1. 

n=0 J \ n=l A 

Then we have for n > 4 

(4.4) /J„(-2Z)+ Y: 2Fu./i,(-2Z) = -e2^°°;i„_4(4Fo), 

k-\-j—n — l 

where Z = (0, Z2, Z3, . . .) with Zn = X! The first equation with n = A gives 

j-\-k—n 

-^Fi3 + 2^^i\ + iF22 = e2^«''. 

This is the first equation of the dispersionless limit of the PfafF lattice equation ( 1 . 18 1 with the 
scaling Fmn 2J^„i„. 



Remark 4.1. The dPfafl curve in (4.2 1 has the structure on (P, Q) e with P ^ and Q ^ q^, 

g2Foo 

which is the same as the dToda curve representing a sphere. The curve has also only two parameters 
of the deformation. However, in terms of the coordinates (A,p) G C^, the curve has an infinite 
number of parameters Fok, unlike the dToda case. In (p, g)-coordinates, the dPfaff curve has an 
asymptotic property 

'^^^^ hi 2mp2™-i ' 

where G2m are determined by G2 and Pgo- Then one can consider the Faber polynomials associated 
with the function q{p). However it seems that there is no role for those polynomials in the dPfaff 
theory. 

5. Universality of the dToda hierarchy 

In this section we first show that the combinatoric descriptions for the functions eo(T) and i?2^''(T) 
in the leading order terms of the free energies (see Theorems 1 1 . 1 1 and 1.2 1 are intimately related. The 



fundamental fact underlying this is that the Euler characteristic two Mobius graphs are equivalent 
to genus ribbon graphs together with some data specifying local orientations at the vertices. The 
relationship between the leading orders of the asymptotic expansions of the free energies log(Z^^) 
in the GUE (/3 = 2), GOE (/3 = 1) and GSE (/3 = 4) random matrix settings, imphes that for 
the special initial condition corresponding to these examples the solutions of the dToda and dPfafT 
hierarchies agree (after a rescaling with proper choice of initial conditions) . We will then show that 
this is merely an example of a more general fact: the solution of the dToda hierarchy is in fact 
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always a solution of the rescaled dPfaff hierarchy, or in other words the dToda hierarchy implies 
the dPfafF hierarchy. However the two hierarchies are not equivalent, the family of solutions to the 
dPfaff hierarchy being much larger (see [221) 

5.1. Universality in random matrix theory. In broad terms universality in random matrix 
theory states that the limiting behavior of eigenvalues of the unitary ensembles of random matrices 
is independent of the external potential Vb(A) in certain scaling regimes. There is also a universality 
for the limiting behavior of the eigenvalues of the random matrices across the orthogonal, unitary, 
and symplectic ensembles in some regions and choices of scalings, see [12l [13] and [30]. In the 
Gaussian case, when Vo(A) is quadratic, we will illustrate universality explicitly by giving a simple 
relationship between the leading order expressions of the free energy F'^''^(ro,T) for the GUE and 
GOE random matrix models (recall that F^^' = F'-'^^ that is, the GOE and GSE cases give the same 
leading order) . A related result is that there is a universality between integrals over the unitary and 
orthogonal groups (not the same expressions as our ziP and Zn"^) after appropriate scalings |ilO, . 

Let us first note the following Lemma relating the leading order combinatoric interpretation given 
by Theorems [TT] and [O} 

Lemma 5.1. Mobius graphs of Euler characteristic 2 are in 1-1 correspondence with pairs made 
up of a ribbon graph of genus together with one of the 2^^^ possible configurations for the local 
orientations of the vertices, where v is the number of vertices of the graph 

Proof. It is easy to see that a ribbon graph of genus together with a choice of local orientations 
of the V vertices gives a Mobius graph of Euler characteristic 2. For a Mobius graph to have Euler 
characteristic 2 there must be a choice of the local orientations of the vertices so that the ribbons 
are untwisted (as the graph lays flat on a sphere), thus giving the pair of a ribbon graph of genus 
and a choice of local orientations. □ 

From this lemma, we see that (for leading order only) there is the following relation between the 



coefficients Ko(j) and /C2(j) of the expansions of the free energies (1.141 and (1.20) 

/C2(j) = r-'^oO), 

where v — "^ji is the number of vertices. We now have: 

i 

Proposition 5.1. The free energies F^^'>{To,T^) have the following relations, 

F(i)(To,T) = fW(To,T) = ^F(2)(2To,2T). 



Proof. The duality between F^^^^ and F^^^ has already been shown in section 1.4 The relation 



between F^-^^ and F^^^ is a direct consequence of the relation between Ko(j) and /C2(j) above. Namely 



from (1.20), we have 



fW(To,T) = E^2(j)(2To)/^ + C^^)(To) 



j! 



E^o{i)2^-\2Toy— + C^'\To) 



J 



J! 



= I E«o(j)(2ro)^'^^ + ^cf (2ro) 

= iF(2)(2ro,2T). 

□ 

This Proposition implies that the free energy of the dToda hierarchy gives an universality over 



those Gaussian ensembles in the leading order with the scalings given in (4.1 1. Also notice that the 



scalings in (4.1 1 explicitly appear in this result. 
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Remark 5.1. It is interesting to note that the scaHngs in Proposition |5.1[ which give universahty 
between the dToda and dPfafF hierarchies, agree with the scahngs found in |IOj giving the universahty 
at leading order between the matrix integrals over the orthogonal and unitary groups (see Theorems 
7.1 and 9.1 of |TQ|)- In our case the partition functions are derived from integrals over the spaces of 
symmetric and hermitian matrices. 

5.2. The dToda hierarchy implies the dPfafF hierarchy. We now show that a solution of the 
dToda hierarchy (2.1) gives a solution of the dPfaff hierarchy (4.2|. This extends the previous results 
for special choices of the free energies associated with the Gaussian ensembles, and it implies that 
the universality holds for the general choice of the potential Vo(A) or the measure in the random 
matrix theory. 



Theorem 5.2. The free energy F of the dToda hierarchy (2.11 also satisfies the dPfaff hierarchy 



Proof. We first calculate the square of the first equation of the dToda hierarchy. 



^2D(\)D(p)F 



2^D(X)D{^L)F 



2 1 



p(A)p(m) p(A)2p(Ai)2 



2D{\)D{^L)F ( ^ _ 



Then from Lemma 2.3 of the kernel formula of p(A), we have 



^2D{\}Dit,}F I ^ 



= 2 



p(a)Mm)' 
p{\)-p{fi) 



- 1. 



Now recall that p{X) for the dToda hierarchy is expressed by p{X) — X — Fqi — D{X)Fi. Then noting 
that 2p{X) = q{X) + A — Foi, we obtain the first equation of the dPfaff hierarchy (|4.2|). 



Now we show that the dPfaff curve q = p 



G2 + e'^^""p ^ is reduced to the dToda curve. If 
are the solutions of the dToda hierarchy, then the G in the curve is given by 

G2 = -4Fn - F^^ + F02 = -2Fn = -2e^°°. 

2Fii and Fn — e^°° (see section 2.1 1. This implies 

o2Foo 



F„ 



Here we have used Fi 



02 



=p^ - 2e 



-Foo 



p^ 



= \P 



p-Fbo 



P 



With g = A — Fqi — 2D{X)Fi — 2p — X + Fqi, this equation leads to the dToda curve, which is 
equivalent to the choice, q = p — ^-y^- □ 

Finally, we remark that this universality can be seen in the dPfaff equation with the rescaling 
and the proper normalization: The dPfaff equation is given by 



2F?, 



which can be split to 



22 



Then if F satisfies the dToda hierarchy, then both sides of this equation are zero. Notice that the 
left hand side is the dKP equation and the other side is the dToda equation. This splitting can be 
also seen in the Pfaff lattice equation (1.15), i.e. the first equation in (1.15) can be written in the 
form, 



_d_ 



-4 



du 



6u 



du 
dTi 



^d^u 



.,d\ 2 



4z;+t; ) . 
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Then one has a KP solution for u which also satisfies the Toda lattice equation given in the right 
hand side (see e.g. [B]). 

6. Appendix A: The dKP hierarchy 



The purpose of this Appendix is to show how one can get the kernel formula (2.3 1 for the KP 
hierarchy, which is one of the most important formulae in this paper. We begin by giving a quick 
introduction to the dKP hierarchy: Let A be a Laurent series in p given by 

A=pH h^H , 

p p^ 

where Uk are the functions of an infinite variables (Ti,T2, . . .). Then the dKP hierarchy is defined 
by the set of equations, 

d\ 

— = A} n=l,2,..., 
where <i>„(A(p)) is the Faber polynomial given by 

<i>„(A(p)) = [A(p)"]>o. 

Here [f{p)]>o is a polynomial part of the Laurent series f{p) in p. In particular, <i>i(A(p)) — p, and 
we denote sometime Ti = X. The bracket {/,<?} is the Poisson bracket defined by 

if \^ 5g _ df dg 
^■''^^ dpdX dXdp' 



The Faber polynomial can be written in the Laurent series (see Definition 2.11, 

<i>„(A) = A" - £ -^Q„,„ . 

m=l ^ 

In particular, we have 

oo 1 

$i(A(p))=p(A) = A"- E T^Qi™- 

m=l 

Then, following the computation in [9J (see also [5]), we first multiply A"~^|^ to p(A), i.e 

P(A)A ^ = A 2^ QimA . 

dp dp „=i dp 

Taking the projection of this equation on the polynomial part. 

Pit = ir — rT*"+i - S Qi^-. 



dp \n J dp \n+ 1 J rn=i dp \n — ni 
Now multiply by fi""' and sum over n > 1, we have 

V — — = V — — f^^\ - V "v ^ ^ 



„=i/i"5pV"/ ri=i M" <9p V'^+ 1/ n=i m=i mfj,"^ ^i'' dp \n - m 



This then gives. 



y ± A (^^^ _ 1 _ V ^ y liL (^^ 



^ A*" y n=i 5p ^ n 

1 a 



Thus we have 



P) Qp y ^ 



1 ^ g / ^n(A(p)) \ 
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Integrating this with respect to p, and fixing the boundary condition at /i = c», we obtain 



exp 



1 



ci>(A) 



exp 
1 

^— ^exp 



m— 1 ^ 
oo oo I/O 

n=l m=l A /i n 



Then expressing the Grunsky coefficients in terms of the free energy, i.e. Cnm — -Qm 
we obtain the dKP hierarchy (2.3 1, 

p{^l)-p{\) ^ ^D{X}D(t,)F 

/i — A 

One should note here that the dispersionless Hirota equation for Fnm is completely determined by 



the function A = p + ^ 



(spectral curve). Then the dKP hierarchy defines an integrable 



deformation of the function A(p), and as a result, the Grunsky coefficients can be expressed by the 
second derivatives of the free energy with the deformation parameters. 

We also remark that in [31] (see also [21]), Takasaki derives this formula as a classical limit of the 
differential Fay identity, which can be obtained from the bilinear identity 

|^e«*'-*-)r(t' r(t + [z-i])=0. 
where £.{t,z) ~ ^ 2"t„ and [z^^] is given by 

n=l 

•f J l_ 

y 2^' 3^' ■ 

Takasaki also derives the dPfaff hierarchy from the classical limit of the differential Fay identity 
obtained from the bilinear identity for the DKP hierarchy [20]. This is remarkable, since there 
seems to be no proper classical limit for the Pfaff lattice hierarchy in the matrix representation 
given in P", ^T]} (this is contrary to the Toda lattice hierarchy) . 

7. Appendix B: Computation of Co(ro) 

Here we give the computation of Cq'^''(To) for /? = 1,2 and 4. We use formula (3.3.10) given in 
[?f] . reproduced here: with 



Zlf'^ ( f A^; ) = / dAi . . . / dK n - A, r exp ( -(^ E A 



one finds the following formula, which we call Mehta's formula. 



Z(« (^^A2;o)=(27r)tr5-'^"^ 

7.1. GUE case. For /3 = 2, we have 

C(2)(To;iV)=log 
From A/fehta's formula, we have 



r 1 



z. 



(2) 



N 



A^;0 



Pi 

2 



Z(2)(A2;0) = (2^)t2-^n^'! 
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A straightforward change of variables gives 



^^^(fA^;o)=iV--(27r)tn,-!. 



Inserting this into our formula for C(To) we find: 



Tl Tl \ 

■ — log iV + - log(27r) + E log(i!) • 



So the problem is to find the leading order behavior of ^ log(j!). We compute this in a generalization 
of the Riemann sums approach to deriving Sterlings approximation. 

We start by rewriting the expression X]log(j!) as a double summation which then simplifies: 

n n j n n 

Eiog(j!)= E Elogfc - E Elogfc 

j=l j—lk—l k—lj—k 

n 

= E ('^ ~ ^ + 1) log k . 

We then approximate this sum by an integral using the trapezoid rule 

E logO'O = 7^^ogn + J{n- x + l)\ogxdx + E , 
j=i ^ 1 

where the error in the approximation can bo computed using the Euler-MacLaurin formula. For our 
purposes it is enough to show that E = O(nlogn) in truth it performs much better. The integral 
then gives 

n 3 1 

E log(i!) = —7^^ + -"-^ logn + ©(nlogn) . 
Inserting this back into our expression for C^^\To; N) we find 

C(')(To;N)=l-n'logn-l ^ ^ 



^n^ log N + -n'^ + 0{n log n) . 



Then in the limit, we have with n = NTq, 



cP (To) := lirn^ ^C{To; N) = ^-T^ log(To) - 1^ ^ 



7.2. GOE case. In the case of /3 = 1, we have 



C«(To;7V)=log 



Mehta's formula then gives 



4nM§A^0)=(27rrr 



7(1) 

^2n 



N 



ri n — 1 



-2n 2n , 

nr(i+i 

i=i ^ 

/tt (2j + 1)! 
22i+i j! 



22"7rtn! []2-2J(2j + l)!. 
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Likewise in this case we find that 



C(i)(To;7V) =log 



Z. 



(1) 

In 



N 



{n^ + ^) log ^ + log 



j(i) 



= n2log-+log(n!)-21og2 ( E J ) + E log[(2i + 1)!] + O(nlogn) . 



Stirlings approximation implies that 

n 

log(n!) = E log A; = O(nlogn) . 

k=l 

So we only need to find the behavior of Elog((2j + 1)!) for large n: as above we follow the basic 
method of Stirlings approximation and write the sum as the approximation of an integral plus an 
error term. For large n the error term is of order 0(1); 

n-l n-12j + l 2n-l n-1 

E log((2j + l)!)= E ElogA:= E E logfc 



j=l fe=l 



Then we find 



which gives 



= E [(n - log(20 + {n-l) log(2l + 1)] 
(=1 

1 n—l 

= -log((2n- 2)(2n- 1)) + / {n - x)log{2x{2x + l))dx + E 

2 1 
3 

= --n^ + n^log(2n) + O(nlogn) . 



C(i) (To; N) = log(2n) - log N - -n^O{n log n) , 
c(^)(To) := ^lirn^ l^C^^)(To; N) = T^\og{2To) - ^T^ . 



7.3. GSE case. For /? = 4, we have 

C7(4)(ro;Ar) = log 
Mehta's formula with /3 = 4 gives 



Z^)(A2;0)=7rt2-2»^+S[](2j)!, 

which together with a change of variables gives 

iV, 



Therefore 

(To; TV) = -n^ log + E log(2j)! + 0{n log n) . 

So we must compute the asymptotic behavior of Elog(2j)! it is not suprising that it is identical 

to the asymptotic behavior of Elog(2j + 1)! found in the /? = 1 case. This is in fact yet more 
evidence of the duality that exists between /3 = 1 and 4. Inserting this back into our expression for 
C^^\To]N) we find 



C('')(To; N) = -n^ logN - -n^ + log(2n) + ©(nlogn) 
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from which wc conclude that 
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